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ABSTRACT 

With the increasing number of directly imaged giant exoplanets the current 
atmosphere models are often not capable of fully explaining the spectra and 
luminosity of the sources. A particularly challenging component of the atmo- 
sphere models is the formation and properties of condensate cloud layers, which 
fundamentally impact the energetics, opacity, and evolution of the planets. 

Here we present a suite of techniques that can be used to estimate the level 
of rotational modulations these planets may show. We propose that the time- 
resolved observations of such periodic photometric and spectroscopic variations 
of extrasolar planets due to their rotation can be used as a powerful tool to probe 
the heterogeneity of their optical surfaces. In this paper we develop simulations to 
explore the capabilities of current and next-generation ground- and space-based 
instruments for this technique. We address and discuss the following questions: 
a) what planet properties can be deduced from the light curve and/or spectra, and 
in particular can we determine rotation periods, spot-coverage, spot colors, spot 
spectra; b) what is the optimal configuration of instrument /wavelength/temporal 
sampling required for these measurements; and, c) can principal component anal- 
ysis be used to invert the light curve and deduce the surface map of the planet. 
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Our simulations describe the expected spectral differences between homoge- 
neous (clear or cloudy) and patchy atmospheres, outline the significance of the 
dominant absorption features of H2O, CH4, and CO and provide a method to 
distinguish these two types of atmospheres. Assuming surfaces with and without 
clouds for most currently imaged planets the current models predict the largest 
variations in the J-band. Simulated photometry from current and future instru- 
ments is used to estimate the level of detectable photometric variations. We 
conclude that future instruments will be able to recover not only the rotation 
periods, cloud cover, cloud colors and spectra but even cloud evolution. We also 
show that a longitudinal map of the planet's atmosphere can be deduced from 
its disk-integrated light curves. 

Subject headings: giant planets: general — direct imaging: surface features, 
photometric and spectroscopic variability, lightcurve inversion 



Introduction 



Clouds play a fundamental but complex role in the atmospheres of brown dwarfs and ex- 
oplanets: describing their vertical and horizontal distributions, composition, formation and 
evolution are the outstanding challenges faced by models of ultracool atmospheres. Refrac- 
tory minerals, such as Ca-Al-oxides, metallic iron and silicates, condense in the temperature 
range ~1, 300-1, 900 K, forming clouds t hat dominate the atmospheres of L- type dwarfs and 
exoplanets of similar temperatures (e.g. iBurrowd l2009l : iFortney et al.ll2008l). These silicate 
clouds are responsible for their very red near-infrared colors (e.g. iBurrowsl l2009t iBurgasser 
2009 : Marlev et al.l 201o[). S ilicate grains, for example, have already been observed with 



Spitzer (ICushing et al.l l2006l ). supporting this assumption. In contrast, the spectrum of 
cooler T-dwarfs is markedly different — it is characterized by blue near-infrared colors and 
is dominated by the absorption of stable gas-phase H2O and CH4. The spectra of T dwarfs 
are explained by clear, i.e. non-cloudy, models. The depletion of refractory elements at 
the L/T-transition regime (a t temperature s of ^1,000 to 1 , 200 K), caused by their "rain 
out" to deeper, hotter layers ( lBurrowsll2009t lBurgasserll2009l : lAUard et al.l 1201 ll ) results in a 
strong change in the spectrum which, combined with the cloud dispersal assumption, is the 
proposed mechanism for the onset of the T-dwarf regime. 

Although the loss of cloud opacity at the L/T transition is qualitatively consistent with 
the drastic changes observed in the spectra, the process leading to the loss of clouds has 
not yet been identified. At least two different ideas have beei i explored. One p l ausibl e 
mechanism, supported by the work of iKnapp et al.l ( 120041 ) and iTsuji fc Nakajimal ( 120031 ). 
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proposes a thinning (or even complete dissipation) of the global cloud cover, caused by the 
growth and subsequent sinking to deeper layers of the cloud particles . Another possibility, 
motivated in part by visible and near-in frared observations of Jupiter (IWestphal et al.lll974j : 



Orton et al. 



19961: 



Dyudina et al.ll200l[ ) and by near-infrared Cassini/VIMS observations of 

Saturn (jBaines et al. 2005h. proposes a sudden appearance of c lear, optically thin "holes" in 



the global cloud deck ( iBurgasser et al.ll2002l : lMarley et al.ll2010l . and others). In this scenario, 
in objects close to the L/T transition flux from deeper, hotter regions can escape through 
these holes. Therefore the atmosphere will appear patchy with bright, hot and deep regions 
next to cooler, optically thick higher-altitude ones. Both scenarios can produce similar 
near-infrared colors at the L/T transition by varying different parameters - temperature, 
surface gravity, sedimentatio n efficiency and/or grairi size for the former and cloud-cover 
and distribution for the latter (jSaumon &: Marleyll2008l ). The predicted differences are subtle 
which, when combined with the uncert ainty in the naeasur ed temperature, currently makes 
it difficult to distinguish these models (IMarley et al.l 120101 ). For example, models that can 
explain the observed broadband photom etry of L/T brown dw arfs do not guarantee a proper 
match for their near- infra red spectra (Burrows et ahlboOGl ). This was clearly shown for 
the case of HR8799b (e.g. Barman et all 201 la ), where the broad band photometry is well 
matched to an L6 dwarf while the near-IR spectrum is not. 

The differences found between the spectra of individual L/T transition objects and the 
model predictions, while difficult to explain, are not surprising as the current models do not 
account for the possibly complex three-dimensional atmospheres and large-scale atmospheric 
circulations, which also set the appearance of Jupiter. Therefore, explaining the atmospheric 
properties of the L/T dwarfs and giant exoplanets of similar temperatures, will likely require 
more than one- dimensional information. Understanding the spectra of these ultra cool atmo- 
spheres will require a physical model not only for the formation, evolution, and destruction 
of cloudy regions, but also for their longitudinal/latitudinal and vertical distributions. 

Not surprisingly, clouds also represent a key problem in the atmospheres of giant ex- 
oplanets of similar temperatures. Many of the recently discovered, directly imaged exo- 
plan ets fall in the temp erature regime of the L/T transit ion: HR8799 bcde wi t h T~1,000K 
e.g. iMarois et al.l |2010[) and (3 Pictoris with T~1,500K ( iLagrange et al.l 120091 : iQuanz et al. 



20101: 



Skemer et a. 



Seiff et al. 



20121) . These planets are cooler than Hot Jupiters (~ 1,500 — 2,400 K, 

Seager fc DemiiiEj 201ol ). but much hotter than the effective temperature of Jupiter (~ 200K 



19981 ). While directly-imaged planets may differ from brown dwarfs in bulk 



chemical composition, surface gravity and formation mechanism, their atmospheric physics 
is thought to be very similar. 



The need for understanding cloud properties became even more pressing with the real- 
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ization that many directly imaged exoplanets are significantly under-luminous in the near- 
infra red bands compared to field brown dwar f s and some sta t e-of-th e -art models: HR879 9b 



(e.g. 



Marois et al.l (120101) 



2007 



Currie et al.l (120111): Barman et al.l fl2011ar): ISkemer et al.l fl2012[)): 



Patience et al.ll2010l : ISkemer et al.ll201ll ). Interestingly, /3 Pic 

appears to alig n 



2M1207b flMohantv et al. 

b - a younger planet with est imated mass close to HR 8799b and 2M1207b 



well with model predictions (ILagrange et al.ll2010l : iQuanz et al.ll2010l : iBonnefoy et al.ll201ll ). 
In contrast, the prominent underluminosity of HR8799b and 2M1207b sparked intense the- 
oretical work and most groups were led t o propose cloud prope rt ies that differ significantly 
from those assumed for field brown dwarfs ( Barman et al. 20 11 al ibi: iMadhusudhan et al. 2011 



Skemer et al.ll2012l ). Several groups found best photometric matches from a combination of 
thick and thin clouds an d therefore argued fo r patchy atmospheres in the HR 8799 planets 



(e.g. iMarois et al.l ( 120081 ): ISkemer et al.l ( 12012! ) ). For a more detailed discussion on the latest 
results we refer the reader to Section 15. 1[ 

While it appears that directly-imaged planets are not a new class of objects but a natural 
continuation of the L-dwarfs sequence (rather than of the T-dwarfs), these developments 
emphasized the need for developing a more realistic model for cloud properties which, in turn, 
demand new, mult i- dimensional data. Photometric variations due to rotating, spectrally 
heter ogenous objects have been proposed as a probe of the c loud p roperties in brown dwarfs 

2001); Burgasser et al. (j2002 )). Similar observations of 



Bailer- Jones fc Mundtl fll999 



(e.g. 

Earth from the EP OXI spacecraft were used to probe land mass and ocean distributions 
(ICowan et al.ll2009l ). While the first searches for varying brown dwarfs produced several 
tentative detections, the past years brought the detection of per i odic, rotational varia tions 
in L/T transition dwarfs JArtigau et al.lbood : Iciarke et al.ll2008[ iRadigan et al.ll2012l ). For 
example, a newly discovered several L/T transition brown dwarf shows a peak-to-peak near - 
infrared variation as large as 27% in the J-band (IRadigan et al.ll2012t lApai et al.l |2012| ). 
With short typical rotation periods (<10 hours) these sources reveal patchy cloud covers on 
brown dwarfs. 

The best example of rotationally-induced photomet ric variability of a giant planet is, 
not surprisingly, in Jupiter. Using IRTF and HST mosaics iGelino &: Marleyl (120001 ) simulated 
the photometric variability of an unresolved Jupiter due to its rotation. At 4.78 /xm Jupiter 
indeed shows a very strong rotational modulation (up to 0.2 magnitudes), detectable at 0.41 
/xm as well on the level of 0.04 magnitudes. The culprit for these variations is the Great Red 
Spot which manifests itself as a large dark patch in the thermal infrared (where the directly 
imaged giant planets are bright). Thus, while not at the temperatures of the current census 
of directly imaged planets, Jupiter can be used as a reasonable starting point for studying 
their expected optical appearance. 
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The need for characterizing cloud properties, combined with the exciting new results on 
brown dwarf variability, and the strong indications of patchy dust clouds in directly imaged 
planets motivate our study to propose variability of directly-imaged giant exoplanets as a 
means for characterizing their cloud covers. We present here a model for the atmospheric 
appearance of a directly-imaged giant planet that can be easily modified to represent any 
scenario: clouds on a clear atmosphere, clear "holes" on a cloudy atmosphere, a global clear 
or cloudy atmosphere with cold and/or hot s pots. Cold spots cou ld represent cloud thick- 



ness/structrure variations and heterogeneity ( iRadigan et al.ll2012l ) which would change the 



brightness of the respective patch so we included both possibilities. The goal of this paper is 
to explore the observational signatures of such heterogenous atmospheres on directly-imaged 
giant planets. The choice of instrument, wavelength, and cadence of the observations are all 
target-dependent and not obvious. New instruments can be optimized for exploiting this 
technique, but this requires an understanding of the variations expected from the rotating 
giant planet targets. We also use the simulated lightcurves to deduce the longitudinal distri- 
bution of the eigencolors using principal component analyses and to recreate the longitudinal 
spot patterns of the input map. 

Here we provide a framework for identifying specific photometric and spectroscopic 
signatures expected of future directly imaged giant planets. The paper is organized as 
follows. In Section [2] we explain the details of our model. Our results are described in 
Section [3l We present the deduced longitudinal map of the planet in Section HI discuss the 
key points in Section [S] and draw our conclusions in Section IS 



2. Model Description 

We constructed a model to predict photometric and spectroscopic variations as a func- 
tion of exoplanet rotation phase, wavelengths of the observations, and instrument/telescope. 
The key properties of the target are its effective temperature and inclination, and the tem- 
perature, spatial and size distribution of the spots. We use the contrast limits for the current 
and future instruments provided in the literature to simulate the relative photometric and 
spectroscopic accuracy. Our model first generates a 2D spot distribution for the exoplanet, 
determines the rotational modulations in the integrated lightcurve and spectrum as a func- 
tion of wavelength, and then simulates the observations with the selected instrument. We 
explore the possibilities of this technique for giant planets to guide future observations and 
identify requirements for various instruments to inform their development. 

We will first review and discuss the cloud model, then describe the spectral libraries we 
use and discuss the simulated observations. 
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2.1. Exoplanet Cloud Models 



We model the spot distribution of the exoplanet with a combination of an ambient/mean 
spectrum, calculated for the effective temperature of the planet, and a set of elliptical spots 
with different temperatures. The spots can have any size, shape, spatial distribution, tem- 
perature, surface gravity, metallicity and covering fraction. In the following, when we speak 
about the features of our models we will refer to "spots", but clouds will be used in the 
context of the astrophysical objects in ultracool atmospheres. 

As a guide for spot size and shape distribution we use Jupiter, our best yet imperfect 
analogue. As will be shown, the general results of this paper do not depend strongly on 
this choice or the specific cloud shape and distribution. Th e cloud pattern of Jupiter i s 
described by latitudinal bands, correlated to zonal circulation ( IVasavada fc Showmanll2005l ). 
The visual appearance of Jupiter is dominated by ammonia clouds located between 250 
mbar and 1 bar (West et al. 2004). Above them are hazes which are observed in the 
near-UV, while deeper atmos pheric layers are probed outside methane bands in the near-IR 
( iBarrado-Izagirre et al.ll2009l ). 



We set up our initial models to resemble the atmosphere of Jupiter, using elliptical spots 
with an aspect-ratio of 1.5, based on the Great Red Spot (GR S) - with dirn ensions of 12,400 
km and 19,800 km it had an aspect ration of 1.59 in 2006 (IRogersI l2008l ). To emphasize 
the effect such a spot can have on the visual appearance of a giant planet, we no te the 
size t he GRS had a 100 years ago when its longitudinal extent was about 45,000 km (jlrwin 
20031 ) and it covered about 3% of the total surface area of Jupiter. The largest attainable 
size of turbulent eddies, defined as the Rhines length, is a function of the atmospheric wind 
speed and the gradient of the Coriolis force ( IShowman et al.ll2010l ) (Section 3.6, Equation 
35), neither of which are constrained for giant planets. We note that the maximum size of 
spots present in planetary atmospheres is not necessary equivalent to the Rhines length, as 
the nature of such spots may be significantly more complicated than simple cyclonic eddies 
with different temperatures. It is not unreasonable to imagine a rather exotic situation 
where the vertical structure of the atmosphere is such that multiple, stratified layers of 
hazes and/or cloud layers with variable thickness alternate in such a way that they do or 
do not obscure deeper/hotter regions in the atmosphere. Nevertheless, the Rhines length 
provides a reasonable initial scale and we adopt it as the parameter setting the size of the 
largest spots. 

For Jupiter, the Rhines scale is on the order of 10% the planet's radius. Atmospheric 
models of Hot Jupiters, h owever, have suggested the presence of very high wind speeds on 
the order of 1-3 km s~^, (IShowman et al.l 120101 ) (Section 3.3, Table 1) and the possibility 
of much larger Rhines scale (comparable to the size of the planet). With a similar radius 
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to Jupiter and assuming comparable rotation rates but faster winds (between Jupiter's and 
those for the Hot Jupiters), the Rhines length for directly-imaged giant planets will be larger 
than that for Jupiter. As a consequence the largest spots may be larger. Thus, to construct 
our atmospheric map we use a single giant spot covering a 5% fraction of the total surface 
area and a set of additional, smaller spots, distributing their semi-major axes as a power-law 
of D{N) = A X 10~°', where D is an array of N semi-major axes (one for each spot), A is 
a scale-factor (described below) and a is the power-law index. The total number of spots, 
the power-law index, the total spot-covering factor (/^ hereafter) and the spot aspect-ratio 
are free parameters in the model. For illustrative purposes, here we use N = 20 spots with 
a = 2.0 and fc = 10%, in such a way that the largest 5 spots correspond to 85% of the 
"patchy" contribution. The choice of parameters is such that a single giant spot dominates 
most of the signal, a few smaller spots produce smaller, albeit still detectable signatures (as 
discussed below) while the rest of the spots are the tail of the power law and too small to 
be detectable. We use the scale-factor A to scale the sizes of the spots such that the total 
area covered by the sum of all spots is the predefined fc, while keeping their relative sizes, 
as determined by the power-law size distribution, unchanged: 



A 



N 



(1) 



where Si is the surface area of the ith spot. 

To calculate what fraction of the surface the different spot types cover during the rota- 
tion, we first randomly distribute the spots on a sphere then project the hemisphere visible for 
discrete rotational phases and finally measure the rotationally modulated fractional coverage 
for each spot. We use orthographic projection (ISnyderl 119871 ) . a technique that represents 
the actual appearance of a distant planet. It does not preserve the size or the shape of the 
surface features but as we are interested in the disk-integrated lightcurve and not in the best 
cartography o f its surface, this transformation is well suited for our purposes. The projection 
is defined as ( 1Snyderlll987l ): 



X = i?cos(0) sin(A — Aq) (2) 
y = R[cos{(f)i) sin(0) — sin(0i) cos(0) cos(A — Aq)] (3) 

where x and y are the cartesian coordinates on the projected 2-dimensional map, A and 
are longitude and latitude on the sphere, (Aq, 0i) are longitude and latitude of the center 
point of the projection and R is the radius of the sphere, which is unity in our model. 
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An example of the two projected hemispheres of the planet with an inclination of 0° is 
shown in Figured], where the different colors correspond to different spots, as described in the 
figure. From this projected map we extract the rotationally modulated fractional coverage 
of each surface element as a function of the planet's period, phase angle and inclination 
as described above. Figure [2] shows the contributions of the largest spots to the visible 
hemisphere as a function of rotational phase. The giant spot covers 21.5% of the projected 
visible hemisphere, while the next five spots by size cover between 2% and 4% each. We note 
that the relationship between the total surface fraction covered fc and the surface fraction 
of each spot is a function of the inclination of the planet and the size and shape of the spots. 

In the next section we describe how these covering fractions are used to combine the 
spectra of all spots present on the hemisphere facing the observer for each rotation phase. 



2.2. Spectra and Spectral Libraries 



To each unique surface el ement ("spot") we a ssign a model spectrum from one of two 



different s pectral libraries - (IBurrows et al.l |2006| ) (BOG hereafter) models or the AMES 
models of (lAUard et al. I l200lh (AOl hereafter). The free parameters of these model libraries 
are temperature, metallicity, log{g) and cloudy or clear atmosphere. We combine the cloudy 
and the clear models from the same library, keeping the surface gravity and the metallicity 
constant for a given object. We note that this step implicitly assumes that the pressure- 
temperature distribution of each column of gas is independen t of that of neighbor ing columns, 
a good first-order assumption, which is not strictly correct (IMarley et al.ll2010l ). 



We explore objects with broad temperature range around the L/T transition (T ~ 700 
K to 1,400 K), representative of the giant self-luminous planets current and next-generation 
facilities are expected to directly image. Throughout the paper we use model spectra for 
solar metallicity and log{g) = 4.5 for simplicity. Model spectra for different temperature 
regimes for AOl (Cond) and BOG (both Clear and Cloudy) are shown in Figure [31 The two 
Clear models are quite similar, the main difference being the more detailed features of AOl. 
We note that the Cloudy BOG are significantly different from the Clear BOG both in the 
strength and in the shape of the spectra in all three filters shown in the figure, a feature that 
will be discussed in more details in the discussion section. 



Next, the spectra of all spots present on the hemisphere facing the observer are weighted 
by their respective surface cover fraction and linearly combined for each rotation phase as: 
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N 

Fx,tot = J2[Fx, X (4) 

i=l 

where Fx^i and fc,i are the flux density and covering factor respectively of each surface 
element on the visible disk of the planet and N is the number of different surfaces on the 
visible disk. The photometry is calculated from the flux density by first normalizing the 
spectrum by the width of the different filters we explore and then integrating. 



2.3. Simulated Observations 



2.3.1. Targets and Instruments 



As a target we assume a star and giant exoplanet resembling HR8799 c (at a distance of 
40 parsec), with a rotation period of 4 hours (planet's apparent brightness: J = 17.65 mag, 
H = 16.93 mag, Ks = 16.33 mag, Tfjf = 1, OOP K aro und an A5V star with J = 5.38 mag, H 
= 5.28 mag and Ks = 5.24 mag, (IMarois et al.ll2008l )). Correspondingly, the planet-to-star 
flux contrast ratio is ^ 
Ks-band. 



1.2 X 10"^ in J-band, ~ 2.2 x 10"^ in H-band and ~ 3.7 x 10" 



m 



We model four different setups representative to the current and next-generation fa- 
cilities shown in Table [1) an 8m-class telescope with Adaptive Optics (8m AO) represent- 
ing VLT/NACO and Keck AO; an 8m-class telescope with an Extreme-AO (8m ExAO) 
representing VLT/SPHERE, Gemini/GPI, LBT/AO; a 30m-class telescope with Extreme 
AO (30m+ ExAO) representing Giant Magellan Telescope (GMT), Thirty Meter Telescope 
(TMT), Extremely Large Telescope (ELT)); and the James Webb Space Telescope (JWST). 
While these examples do n ot include all p lanned instruments, s uch as the Lyot 1800 project 
(jOppenheimer et al.l |2004| ) or ATLAST (jPostman et al.l l2010l ). they bracket the range of 
relevant instrument capabilities for the next two decades. 

To estimate the achievable photometric accuracy of the different instruments we rely on 
the residual radial contras t curves provided by the inst rument teanis: VLT/NACO at 4 nm 
jKasper et al.ll2007l . I2OO9I) , VLT/SPHERE in J- band Jvigan et allboiol : iMesa et al.l I20 111) 



Gem ini/GPI in H-band jMacintosh et al.l I2OO8I), TMT/ PFI in H-band (Macintosh et Z 



20061). ELT /EPICS in J-band flKasper et al.ll2008l . l2010f ) and JWST/NIRCAM in K-band 
(IGreen et al.ll2005l ). All contrast limits are for coronographic images. For VLT/NACO, 
VLT/SPHERE, Gemini/GPI and JWST/NIRCAM we place the planet at a separation of 



1", similar to the "c" 



planet in the HR8799bcde system (IMarois et al.ll2008l ). For the 30m 
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class telescopes we assume a separation of 0.2", due to the smaller field-of-view of some of 
these instruments. The sensitivity limits for the four different instrument classes we explore 
are shown in Table [2J 



2.3.2. Synthetic Photometry and Spectroscopy 



To explore our ability to recover surface details of giant exoplanets we create a set of 
simulated observations. Rotationally-modulated lightcurves in J, H and Ks-bands produced 
as described in Section for an exoplanet with an atmospheric map from Figured] are used 
to simulate the flux from the planet measured by a suite of current and future instruments. 
Throughout this work we use the VLT/NACO filters defined with central wavelengths and 
widths (in /im) as follows: J-band (1.265 and 0.25 respectively); H-band (1.66 and 0.33); Ks- 
band (2.18 and 0.35), L' (3.8 and 0.62) and M' (4.78 and 0.59). The actual filter transmission 
curves are available on the ESO instrument website. 

Most state-of-the-art ground- and space-based high-contrast observations rely on rela- 
tive instrument-sky rotations to separate faint point sources from instrument speckles. The 
ground-based version of this technique is often referred to as angular differential imaging 



(API ) - a technique with a variety of implemen tations (e.g. iMarois et al.ll2006al : lApai et al. 



20081 : iKasper et al.l 120071 : iLafreniere et al.l 120071 ). While a powerful method, ADI-type ob- 



servations necessarily pose an important constraint on time-series observations. Because the 
observations rely on field rotation to separate real sources from speckles, for each indepen- 
dent photometric measurement a minimum field rotation rate per image (an upper limit on 
cadence) is required. We describe this by requiring the arc traced by the planet during the 
rotation to be larger than at least four times the full width at half maximum of the point 
spread function (PSF) to ensure that the planet's PSF is well separated from instrument 
speckles. To c alculate t he fie ld rotation rate and time required by the ADI we follow the 
prescription of iMcLeanI (Il997l ) (Chapter 3, equation 3.14) as follows: 



cos{A) cos{(f)) 



rot 



sm U 



Sep 



(5) 
(6) 



where u is the field rotation rate in radians per second, Q = 7.2925 x 10"^ radians/sec 
is the sidereal rate, A is the target's azimuth, is the latitude of the observatory, z is 
the target's zenith distance, frot is the minimum field rotation required for ADI in radians. 
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is the size of the PSF in arcseconds and r^ep is the radial distance from the axis of 
rotation in arcsecond. It is important to note that for a 30-m class telescope the gain in 
cadence is proportional to the primary mirror's diameter (see Eq. [6]). This is because the 
diffraction-limited PSF is smaller, which in turns makes the necessary arc traced during 
the ADI rotation shorter. Space telescopes that allow the rotation of the entire spacecraft, 
such as the Hubble Space Telescope, may in theory allow for a better temporal sampling, 
but in practice such rotations are often time-consuming and may limit the cadence. The 
factor of 4 in the nominator of Equation |6] comes from the requirement for non-overlapping 
PSFs. Throughout this paper we use a target with a declination of 5 ~ +21° (reminiscent of 
HR8799) observed from the latitude of the Paranal observatory (0 = —24°). Using Equation 
|5]for these set of parameters we obtain a value for u near meridian of ~ 0.02° /sec 

Due to the complexity of adaptive optics instruments and the sensitivity of the AO 
correction to the atmospheric conditions it can be challenging to reach absolute photometry 
with 1% accuracy. While extreme AO systems are expected to reach very high signal to 
noise ratios on bright giant planets, as discussed in Section 13. 3[ here we also offer three 
further techniques that will help reaching high-precision photometry. First, we point out 
that rotational mapping does not necessarily require absolute photometry. In some cases the 
host star itself can be used as a photometric reference source. This is not always possible 
due to the high contrast but other planets in the system should provide ideal comparison 
points. A good example is the HR8799bcde system, where relative photometry between the 



three , similarly bright outer planets can provide accurate relative measurements ( lApai et al. 



20121 ). Second, most AO high-contrast imaging pipelines allow the injection of an artificial 
star in the raw data, which can quantify flux losses during the data reduction. Such artficial 
star tests can be combined with the relative photometry to further quantify losses and 
uncertainties. Third, we point out that high-ord er deformable mirror s can be used to inject 



an artificial "speckle" into the optical system ( iMarois et al.l l2006bl ). Such "speckles" are 
images of the star and thus can serve as ideal references for relative photometry, even if 
there are no suitable planets or if the star itself is too bright. Therefore, reaching even sub- 
percent accuracy in relative photometry with next-generation AO systems seems plausible. 

We explore six distinct realizations of the different surface types present on the atmo- 
sphere of the giant planet, shown in Table E) a) cloudy spots on a clear surface, with the 
same (Model Al hereafter) or with different temperatures (Model A2 hereafter), representing 
clouds on a cloud-free surface; b) clear spots on a cloudy surface, with the same (Model Bl 
hereafter) or with different temperatures (Model B2 hereafter), representing clear, deeper 
holes in an otherwise global cloud cover; c) clear surface with cold and hot clear spots (Model 
C hereafter); and d) cloudy surface with cold and hot cloudy spots (Model D hereafter). Us- 
ing the B06 cloudy and clear models we simulate spectral modulations due to rotation for 
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effective temperatures between 800K and 1200K and a wavelength range between 1 and 12 
/im. 

In the next section we present the results of our model using a spectral resolution of 
100 and the simple spot distribution in Figure [1] 



3. Results 

The model described in the previous section allows us to predict rotational variations 
in the photometry and spectroscopy of giant exoplanets and simulate their observations 
with existing and next-generation facilities. Here we will evaluate the capabilities of differ- 
ent telescope/instrument classes for characterizing giant exoplanets beyond one-dimensional 
measurements. We also identify the ideal instrument setup as a function of target tempera- 
tures and cloud properties. 



3.1. Properties of the Simulated Variability 

An example of normalized model lightcurves is shown in Figure H] where we plot the 
results for the J, H, Ks and L' filters for Model B2, as described in Section 12.3.21 The 
maximum normalized amplitude in all three bands occurs when the giant, Clear hot spot 
(which is 200 K hotter than the 1,000 K effective temperature of the Cloudy surface) rotates 
into view. As expected from a closer inspection of the green and yellow lines in the top panel 
of Figure |3l the largest photometric variations, up to 19% from the m ean occur in the J- 



band. This behavior is consistent with results from lArtigau et al.l ( 120091 ) and iRadigan et al. 



(120121 ) who also reported the largest photometric modulations in the J band for two early 
T-type brown dwarfs. The H and K-bands have similar behavior in Figure H] with 14% 
maximum photometric amplitude for the former and 15% for the latter while the amplitude 
of the modulation in L' is not more than 10%. The minimum in all bands is caused by the 
two Clear (red spectrum, top panel of Figure [3]) cold (T=800 K) spots seen on the right 
panel of Figure [TJ This lightcurve minimum is prominent - up to 13% lower than the mean 
J band lightcurve - given that neither of the cold spots cover more than 4% of the visible 
hemisphere (see Section 12.11 and Figure [2]) . Note that the trough is not a symmetric feature 
because the giant spot rotates into view and dominates the lightcurve, decreasing the effect 
of the two cold spots. There are many small spots in the model (less than 1% of the visible 
hemisphere covered) that contribute small-level photometric variations, whose presence can 
be deduced from the lightcurve with sufficient sampling rate and accuracy. 
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Simulated normalized rotational spectral modulations for Model B2 are shown in Figure 
[5] for the JWST/NIRSPEC and JWST/MIRI wavelength range with a spectral resolution 
of 100. The variations are very strong (up to 35% minimum to maximum) not only in 
the J, H and Ks-bands but also in the mid-infrared (MIR) where they can be as high as 
20% at 6 ^m. The rotation phases for the rise and fall of the amplitudes in the MIR show 
distinct variations, in particularly at 9.7 fim compared to 10.3 /xm. The heterogeneous 
nature of the atmosphere is easily discernible as spots of different temperature and chemical 
composition rotate in and out of view. Such spectral maps will not only suggest variations in 
the abundance of gas-phase absorbers, in particular methane or water, but will also further 
constrain the covering fraction and longitudinal distribution of the spots. 

To further explore the behavior of spectral maps and to study how they respond to 
changes in the input parameters (surface type, temperature, covering fraction), in Figure E] 
we present three different realizations of the spectral map in Figure [5] for wavelength range 
between 1 and 2.5 /xm. Here we change only the temperature of the spots and their type 
(Cloudy or Clear) but keep the Te// = 1,000K and total spot-covering fraction fc constant at 
10% and using the atmospheric map in Figured] throughout. Changing the covering fraction 
will change the amplitude of the variations, but not the wavelengths/spectra or the periods. 
Using the nomenclature of Table [3l we show Model Al, Model A2 and Model C on the top, 
middle and lower panels respectively. The relative variations in the J, H, and Ks-bands are 
different between the different scenarios. When the giant spot rotates into view, the J-band 
modulations for Model Al are negative (as there is less of the brighter. Clear surface; green 
line in top panel of Figure [3]), while those in H- and Ks-bands are positive (as there is more of 
the brighter. Cloudy surface; yellow line. Figure [3]). On the contrary, the J-band variations 
are positive for both Model A2 (middle panel) and Model C (lower panel) when the giant 
spot faces the observer. The alternating darker and brighter stripes in the J band in the 
middle panel of Figure [H] are a direct consequence of the differences between the 1,000 K 
Clear (top panel Figure [31 green) effective surface and the 1,200 K Cloudy (top panel Figure 
[31 light orange) model spectra - the former is higher but narrower compared to the latter 
in this waveband. When the Cloudy Spot rotates into view, it enhances the flux density 
(white, J-band in the middle panel of Figure ^ in the sides and decreases it in the middle 
of the filter band, where the contribution of the brighter. Clear surface is lower. Therefore, 
while photometric variations like the ones in Figure [Dean be interpreted as caused by either a 
hotter/colder region and/or different chemical composition, sjjectra/ modulations, such as the 
ones seen in the J-band in middle panel of Figure [6l can identify gas-phase compositional 
variations. The magnitudes of the variations are different between the three models: the 
variations are strongest in Model C (lower panel, minimum to maximum amplitude of 20%) 
and weakest in Model Al (top panel, minimum to maximum amplitude of 10%). We note 
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that the magnitudes of the variations for all three models in Figure |6] are significantly smaller 
than those in Figure [5] where the minimum to maximum amplitude is up to 35% (for model 
B2 from l2.3!^ . Additionally, detection of such spectral modulations can also further improve 
light curve inversion techniques such as the principal component analysis discussed in Section 

SI 

Next we explore the optimal filter for detecting variations for Model B2 as a function 
of the temperature of the giant spot. In Figure [7| we show the results for three different 
effective temperatures of the planet - 800 K, 1,000 K and 1,200 K (top, middle and lower 
panels, respectively) for the J, H and Ks filters. The temperature of the giant spot varies 
between — 300K and +300K from the effective temperature with 5T = lOOi^', while the 
temperatures of the rest of the spots (either cooler or hotter by 200 K compared to the ef- 
fective temperature) stay the same throughout. Changing the temperature of the giant spot 
while keeping everything else constant will not only change the amplitude of the variations, 
but also the shape of the lightcurve with the result that different features will produce the 
largest amplitudes. The two cold spots in the right panel of Figure [1] are very prominent in 
some cases where the giant spot plays only a secondary role (for example, Tg/j = 1200K and 
Tspot = 1000, 1100, 1300K) and are in fact responsible for the largest amplitudes. Regardless 
of the features responsible, the largest photometric modulations - the variations producing 
the strongest observable signal - occur in the Ks-band for Tg// = SOOK (up to 60% for 
the hottest spot, upper panel. Figure [7]) and in the J-band for T^fj = lOOOi^' and 1,2007^ 
(up to 30%, middle and lower panels. Figure [7]). This result sugg ests that for sources with 



physics /c hemistry similar to those in Model B2, like SIMP0136 (lArtigau et al.l l2009l ) and 



2M2139 (IRadigan et al.ll2012l ). the most appropriate filter where the photometric variations 
are largest and will be most easily detected is indeed J, as these authors noted. For colder 
sources or sources with a cold spot, however, Ks would be more appropriate. As expected, 
all three bands "brighten" as the temperature of the giant spot incre ases from Tpff + lOOK 



to Tpf f + SOOK, reminisc ent of the well-known J-brightening effect (IBurgasser et al.l 12000 



2002 



Leggett et al.ll2000l ) probably caused by the appearance of cloud-free regions at the 
L/T transition. Here, the brightness increases not due to increase in the size of the clear 
regions, but due to the increase in their temperature. The three bands brighten at different 
rates, from the slowest increase in amplitude of 12% (Ks-band, lower panel) to the fastest 
of 40% (Ks-band, upper panel). The fastest brightening for all three filters consistently 
occurs for the model in the upper panel. We also note the continuous transition in the most 
appropriate choice of filter where the largest amplitudes occur. As the Teff decreases from 
1200K to SOOK the best filter changes from J (lower panel) to Ks (upper panel) and as the 
temperature of the giant spot decreases from T^jf + 300K to T^jf — 300K, the best filter 
changes from J to Ks (middle and lower panels). With even larger temperature differences. 
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Ks-band amplitudes would eventually exceed those seen in the J-band. Such large temper- 
atu re differences betwee n the different surface types are, however, not expected (although 



see 



Skemer et al.l ( 120121 ) ). All the results discussed above are the direct consequence of the 
interplay between the spectral differences between the Clear and Cloudy model for these 
temperatures (see Figure [3]) and the distribution and sizes of the different surface features. 
The largest photometric amplitudes occur in different filters for different models. For ex- 
ample, for Model A2 the largest modulation for the same calculation as done for Figure [7] 
consistently occur in the Ks-band. A survey over targets with different T^fj can populate 
Figure [7] with more points and in combination with the model presented here can provide 
information on this level of details. 

The following four subsections describe the outcome of simulated photometry for dif- 
ferent instruments classes using the theoretical lightcurves from Figure S] and outline the 
increasingly complex picture of the optical appearance of the giant exoplanet that can be 
learned. 



3.2. 8-m Class Telescopes with AO 



The general prediction of our model is that the largest amplitude changes occur in J- 
band. However, J-band is not the optimal wavelength for current 8 m-class telescopes with 
state-of-the-art AO systems: the planet-to-star flux ratio is ~ 1.2 x 10^^ (see Section l2.3.ip . 
Instead, longer wavelengths provide a better AO correction and more favorable planet-to-star 
contrasts. We use he re the example of VLT/NACO and assurne residual stellar PSF in Ks 
similar to that in L' (IKasper et al.l 120071 120091 : iLagrange et al.ll2010l : iBonnefoy et al.l 120111 ) 
to simulate Ks-band photometry where the planet-to-star flux ratio is higher (the planet 
is brighter): ~ 3.7 x 10"^ (Section [2X1]). For the Ks-band VLT/NACO PSF FWHM of ~ 
0.069 " we calculate frot ~ 16° from Equation |6l Using this value for frot and a field rotation 
rate near meridian u ~ 0.02^ results in ~14 minutes minimum required time for proper 
ADI reduction. The signal-to-noise ratio (SNR) for the planet is ~ 15; here we adopted a 
read noise of 46.20 electrons as giv en by the VLT/NA CO Exposure Time Calculator (ETC) 
and an observing mode similar to (IKasper et al.ll2009l ) but modified for Ks, namely dithered 
individual 1 min exposures with detector integration time (DIT) of 0.3454 sec (given by 
the VLT/NACO manual) and 200 DIT per position. The simulated Ks-band lightcurve for 
model B2 with added normally distributed photometric uncertainty of 10% is shown in the 
top panel of Figure El 



Detection of the photometric modulation is not background-limited but is instead lim- 
ited by the residual stellar point spread function. The smallest detectable Ks-band vari- 
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ability at 3cr level is ~ 20%, on par with our simulated maximum amplitude (see Figure 
HJ. Therefore, this instrument configuration is capable of confidently detecting the Ks-band 
photometric variations produced by the giant spot in Figure [T] and retrieving the rotation 
period of the planet. It is, however, inadequate to explore the smaller variations (at the 
~1% level) expected from the smaller clouds. 



3.3. 8— m Class Telescopes with next-generation Extreme AO 

The next-generation of high-order AO systems (extreme AO, or ExAO) will be capable 
of delivering very high quality correction also at shorter wavelengths, allowing variability 
searches to focus on the most favorable wavelengths. To evaluate this configuration we 
repeat our calculation for the minimum ADI cadence but for J-band PSF - in this case frot 
~ 9° and the minimum required time for proper ADI reduction is ~8 minutes. 



Following the science r equirements for VLT/SPHERE and Gemini/GPI (IMesa et al. 



201 ll : [Macintosh et al.ll2008l ) we use a residual stellar PSF of 10^'' at a separation of 1", and 
a read-noise of 10 electrons to simulate the ~ 100 SNR J-band photometry (with added 
normally distributed noise of 1%, see below) for model B2 shown in the middle panel of 
Figure El We note that the high Strehl ratio provided by the Extreme AO (up to 80%) 
reduces the speckle noise to such low levels that the limiting factors for the detection of 
photometric variability will be the photometric precision, read noise and instrum ent stability. 



Here we use a 1% J-band photometric precision following the prescription of IVigan et al. 



(l2010l ) for VLT/SPHERE, thus setting the limits on the smallest detectable variability. This 
instrument configuration will be able to easily measure the rotational periods of planets, will 
allow preliminary exploration of cloud colors and can even detect the temporal evolution of 
the cloud cover. Evidence for such change s in cloud cover in brown dwarfs have been reported 



( lArtigau et al.ll2009l : iRadigan et al.ll2012l ). but their detection requires high signal-to-noise. 



For example, the amplitude difference in the simulated (H— Ks) color for the giant spot 
in the right panel of Figure [T] is on the order of 0.5% while the (J— Ks) amplitude color 
difference is up to 4%. Amplitude difference in the (J— Ks) color for the two cold spots in 
the right panel of Figured] is ~2% while that in (H— Ks) is ~0.5%. Here we expect a few a 
detection of (J— Ks) modulations caused by the rotation of the giant spot but in general the 
confidence with which such measurements could be obtained will be ultimately dependent 
on the photometric precision. A limitation on the capabilities of this class of instruments 
will be the achievable cadence due to the ADI requirements as discussed in Section I2.3.2[ an 
obstacle that will be overcome by the class of instruments presented in the next section. 
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3.4. 30— m Class Telescopes with next-generation Extreme AO (30m+ ExAO) 



This group of instruments represents the suite of planned extremely-large-aperture 
ground-based telescopes. Following the method described in the previous two sections, 
we calculate a J-band ADI rotation time requirement of ~ 2 minutes and ~ 10 minutes for a 
planet at a radial separation from the parent star of 1" and 0.2" respectively; here we p resent 
the results for the latter. Using the required ch aracteristics of the GMT/TMT/ELT (IGMT 
20061 : [Macintosh et al.ll2006l : iKasper et al.ll2010l ). we use a residual stellar PSF of 10~* at the 
radial separation of 0.2" to simulate the very high SNR (on the order of several thousands, 
using the ELT ETC provided by ESO) J-band lightcurve for model B2 in the lower panel of 
Figure m As was the case for the 8m ExAO systems, the residual stellar PSF is no longer the 
main culprit for systematic errors. H ere we adopt a photorn etric precision of 1% following 



the discussion in iDekany et al.l ( 120041 ) and iLiske et al.l ( 120111 ) . 



The capabilities of these instruments represent a significant step forward from those of 
the current 8- and 10-m class telescopes. The high Strehl ratio combined with the improved 
cadence provided by the 30-m class ground-based telescopes will allow precise measurements 
of the basic characteristics of the planet's atmosphere like rotation rate or large-scale in- 
homogeneity. As for the 8-m class Extreme AO instruments, the 30-m class Extreme AO 
instruments will also be able to measure cloud colors only down to a level ultimately lim- 
ited by systematic trends in the photometry, due to the instrumental changes, atmospheric 
variations, etc. The most significant improvement of the 30-m Extreme AO over the 8-m 
Extreme AO instruments will be in their ability to do moderate-resolution (R ~ 100 - 1000) 
spectroscopy both in the near IR with an Integral Field Spectrograph (IFS) on ELT /EPICS 



TMT/NFIR AOS and GMT/NIRExAO Imager (IVerinaud et all l2010l : iHerriot et al.l 12006 



Johns! |2008[ ) and in the mid-IR using, for examp le, instruments like GMT/TIGER in the 
8 to 20 /zm wavelength coverage (jJaffe et al.l 120101 ). Further improvements will come in the 
form of polarimetric measurements and the extended time-domain probed with the faster 
cadence, allowing detection of smaller-scale variations. Also, the radial separations at which 
giant exoplanets will be observed with a contrast of up to 10^^ will be as small as 0.1", 
allowing the observations of colder and /or older giant pla nets in reflected light at orbital 
distances comparable to that of Jupiter (IKasper et al.ll2010l ) and of very young planets close 
to the snow line. Access to such small angular separations will also significantly expand 
the available distance over which systems containing self-luminous Jovian planets can be 
directly imaged. This illustrates the unique capabilities of the next-generation extremely- 
large telescopes as even the future space-based telescopes like the JWST (discussed in the 
next section) will not have such small angular resolution. Direct imagin g of Neptune-size o r 
even Super Earths around the closest stars should be also achievable ( iKasper et al.ll2010l ). 
The favorable contrast between such planets and their host star at the mid-IR wavelengths 
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will, albeit through challenging observations, allow low-resolution spectral characterization 
of their atmospheres by detecting H2O, CH4, CO, CO2 and NH3 features. Such measure- 
ments could also discern between a hydrogen-rich and a water-rich atmosphere and test for 
presence of hazes and non-equilibriurn cheni is try in the atmosp h eres o f directly-imaged pla n- 
ets (e.g. iMiller-Ricci Kempton et~aD fcoilh : iRoeers fc Seageil fcoioh : IPesert et alTfcoilh l 
Finding such planets in the habitable zone and studying their spectra for biosignatures will 
also significantly advance the field of astrobiology. 



3.5. JWST 

The James Webb Space Telescope will open a new chapter in the studies of directly- 
imaged giant exoplanets. With an aperture and Strehl ratio on par with the VLT/SPHERE 
but free of the atmospheric aberrations, and combined with a dramatically low thermal 
background, JWST will have a greatly increased sensitivity and cadence. To simulate J- 
band photometry f rom NIRCAM we use a residual stellar PSF of 10^^ at a separation of 1" 



using the result s of iGreen et a 



deconvolution ([Gardner et al 



■ (I2OO5I) . To reduce the effect of aberrations, we assume a roU- 
20061 ) mode to be employed by the JWST using a conservative 



caden ce of 15 minutes. Also, following the discussion in iLagage &: European MIRI Team 
(I2OIOI ). we assume a photometric precision of ~ lO"'^. This is already at least an order of 



magnitude higher than the one for the ground-based facilities discusses in Sections 13.31 and 
I3.4l but is still significantly higher than the achievable contrast, and is again what ultimately 
limits the level of detectable variations. 

The simulated observation is shown in Figure M The significant decrease in the speckle 
pattern, combined with the very low background noise, the stability of the instrument and 
the short cadence will allow detailed studies of the rotation periods, cloud distribution, cloud 
colors and spectra, atmospheric maps and possibly weather patterns of directly imaged giant 
planets. The unique strength of JWST will be in the complementary observations in the 
near- and mid-infrared, studying wavelengths that are difficult to explore from the ground, 
such as the pea k in thermal emission of young Jovian planets at 4.5 /im. As discussed by 



Clampiru (120081 ) , both the capabilities of NIRCAM in the F460M filter at separations larger 
than 0.6" and the sensitivity of MIRI will exceed that of even the planned 30-m class ground- 
based instruments. MIRI will also provide a platform uniquely suitable for observing planets 
at wavelengths longer than 5 /xm. 

The fact that JWST will not suffer from telluric absorption will also allow more accurate 
observations of several key gas-phase molecular absorbers. As seen in Figures [5] and [6] 
prominent absorbers such as H2O, CH4, CO and NH3 with very strong spectral features can 
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be studied by NIRSPEC and MIRI in detail, allowing their spatial and tem pora l distribution 
among other things to be deduced. As discussed by iMarley et al.l (I2OIOI ) and iBailer- Jones 
( 120081 ) detection of variability in very temperature-dependent spectral features such as FeH, 
Nal and KI and/or correlated modulations between these features will further support the 
assumption of a patchy atmosphere. Detection of strong water absorption, for example, will 
indicate a hydrogen-rich envelope. The JWST will also be uniquely suited to study the 
mid-infrared regime, where the spectroscopic variations can be as high as 20% for the model 
seen in Figure O 



4. Lightcurve Inversion 

To explore the feasibility of recovering the map from Figure [H we use a principal com- 
ponent analysis (PCA) and sinusoidal variations in surface brightness as a function of only 
the longitude as studied by Cowan et al. (2009) in combination with the formalism for light 



curve inversion (LCI) developed in iRusselll ( 1l906l ). The purpose of this section is to illustrate 
the potential of LCI. 

As an example we utilize a covariance-matrix PCA, using the five J, H, Ks, L' and M' 
filters as the input vectors defining a 5-dimensional parameter space to be converted to a 
minimum set of eigenvectors sufficient to produce the simulated variations. The resulting 
normalized eigenvalues are shown in Figure [TOl Following the prescription for the method, we 
notice that the largest gap occurs between the second and the third eigenvalues, suggesting 
that only two principal colors above the mean background color are required to reproduce 
the observed variations, consistent with the three colors we used for our map in Figure [TJ 
The contribution from the two primary colors as a function of the rotation phase is shown in 
Figure [11] where we use two panels to show the actual magnitude of the variations (rotation 
phase of zero is defined as the left panel on[T]). A maximum value corresponds to a rotation 
phase at which the largest area covered of that respective eigencolor is present on the side 
facing us and a minimum - the least. The primary eigencolor basically traces the input 
(orange) J-band lightcurve from Figure HI On the contrary, the behavi or of the secondary i s 



not immediately obvious from the input lightcurve. As pointed out in ICowan et al.l ( l2009l ). 
the eigencolors do not represent an actual color but a deviation from the mean "color", which 
for our map in Figured] is the 1,000 K (orange) background temperature. 



To recover the surface map from Figured] we follow the prescription of iRusselll ( 119061 ). 
We use the two different eigencolors as separate light curves and expand them in spherical 
harmonics. We then invert them for a longitudinal, sinusoidal brightness distribution. The 
inversion is done in one dimension only - there is no latitudinal resolution. For simplicity. 
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we assume that the planet's rotation axis is perpendicular to the line of sight. The resulting 
longitudinal maps in Mollweide projection for the primary and secondary eigencolors are 
shown in the lower panels Figures [T2] and [T31 respectively. The " other" eigencolor on each 
map has been " masked out" because the two eigencolors are orthogonal. In the upper panels 
of Figures [12] and [13] we compare the PCA-inverted longitudinal maps to the input map 
(in the same projection) from Figure [1] The two figures show that this simple method 
successfully recovers the longitudinal position of the major features, namely the giant hot 
spot at a longitude of +135° (Figure [T^]) and the two groups of smaller, cooler spots at 
longitudes of —60° and +60° respectively (Figure [T3|) . 



As iRusselll (119061 ) pointed out expanding the lightcurve in spherical harmonics does 
not completely determine the harmonics expansion of the surface brightness, as all but the 
first odd harmonics are absent in the lightcurve. Nevertheless, this simple model not only 
recovers the contribution to the lightcurve from the primary, dominant eigencolor (maximum 
variation of up to 0.3, upper panel. Figure [TT|) but also identifies the effect of the much 
weaker components (maximum variation of 0.013, lower panel. Figure [TT]) from the secondary 
eigencolor, allowing us to infer the presence of distinct surface features on the planet. 

A more comprehensive way to obtain the surface map of the plane would be to use the 
PCA on the spectroscopic instead of the photometric variations, using the resolution of the 
instrument to increase the dimensionality of the data - a single, low-resolution spectrum of 
R~ 50 will increase the number of orthogonal data sets by an order of magnitude over the 
five filters discussed above. 

The method described thus far is well suited to study the color asymmetry caused by 
the spots in the form of spatial distribution over the surface of the planet, aft er simplifying 



assumptions like the longitudinal sinusoidal intensity map or an N-slice map (see I Cowan et al. 



20091). 



5. Discussion 



5.1. Patchy Dusty Atmospheres in Brown Dwarfs and Exoplanets 



Early attempts to understand the nature of giant planets were based on the use of 
atmospheric models of L and T brown dwarfs as templates to assess the detection sensi- 
tivity for direct imaging surveys. While at similar temperatures to the regime of the L/T 
transition objects (between ~ 900K and 1200K where dust clouds settle below the photo- 
sphere), it has been sugge sted that giant planets may, in fact, remai i i cloudy where the T 
dwarfs become cloud-free ( ISaumon fc Marleyll2008l : ICurrie et al.ll201ll : [Barman et al.ll2011al : 
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Madhusudhan et al.ll2011t ISkemer et al.ll2012l ). reminiscent of the results of [Stephens et al. 
( 120091 ) on low-gravity BDs. This does not necessary invoke a new kind of atmospheres for 
the directly-imaged giant planets. Most likely it is an age effect of atmospheres rather more 
dusty for their spectral type. 

The observed abundance of methane compared to CO in the four planets of the HR8799 
system has been reported to be lower compared to chemical equilibrium models predictions, 
sugg esting efficient convective mixing on time-scales shorter tha n the CO to CH4 equilibrium 
rate (IHinz et al.ll2010l : iBarman et al.ll2011al : ISkemer et al.ll2012l ). This assumption is further 
supported by the 3.88 to 4.10 utsi spectroscopy of HR8799c being inc onsistent with both 
the DUSTY dchabrier et"aDl200oh and the COND jBaraffe et al.lboosh models, indicating 
non-equilibrium chemistry at work. It appears that giant planets such as HR8799bcde and 
2M1207b (a rather more pronounce d case) may indeed have the temperature of T-dwarfs but 



the du sty appearance of L-dwarfs (jChauvin et al.l l2004l : ISkemer et al.l I2OIII : iBarman et al. 



2011bl ). In particular, the luminosity of HR8799b, for example, implies an effe ctive temper 



ature of 850K while its colors suggest a much hotter atmosphere of 1300K (ISkemer et al. 



2OI2I ). Models with thick clouds (with the cloud base sig nificantly dee p er com pared to that 
of L/T dwarfs) and low surface gravity are invoked by ICurrie et al.l (120111 ) to match the 
planets in the HR8799 system. Their best-fit model for HR8799b yields T^fj = 800-lOOOK 
and log(g) = 4.5 and Te// = 1000-1200K and log(g) = 4-4.5 for HR8799cd. The authors 
put less importance on the effects of non-equilibrium chemistry in reproducing the 1 to 5 
Hm SEDs but suggest that the planet atmospheres may indeed be out of equilibrium. Inter- 
estingly, they also argue that "patchy" cloudy models may provide an even better fit to the 
data. On the contrary, to explain the observed low effective temperature (less than lOOOK, 
ty pical of cloud- f ree T-d warfs), red colors and smooth spectrum of HR8799b, the results 
of IBarman et al.l (l2011al ) support the presence of thick photospheric clouds and enhanced 
metallicity in the presence of non-equilibrium chemistry. The authors note that the spectrum 
of the planet is markedly different from that of typical field brown dwarfs, namely in the 
weaker methane and CO and strong er water absorption (suggesting a hydrogen-rich atmo- 
sphere). In contrast to the results of lCurrie et al.l (201ll) who st ress the importance of cloud 
thickness over non-equilibrium chemistry, IBarman et al.l (l2011al ) assign equal importance to 
both. However, their best-fit model for solar abundance, with Tg// = HOOK and log(g) = 3.5, 
indicates a radius that is smaller than theoretically expected, suggesting that the metallicity 
of the planet must be higher. A chemical enrichment by a factor of ten produces a signifi- 
cantly better fit for their data to a model with Teff = 869K and log(g) = 4.3, consistent with 
the planet's mass and age. Similar scenario can be drawn for the case of 2M1207b, where 
the apparent low luminosity and red colors are in stark contrast with the spectrum-derived 
Te//~1600K. Using an atmospheric model with thick clouds, non-equilibrium chemistry, a 
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mixture of grains and low surface gravity, iBarman et al.l ( 1201 lb[ ) propose a best-fit solution 



with T 



lOOOK and log(g) = 4., consistent with cooling track predictions and disfavor- 
ing exoti c scenarios such as ed ge-on disks or planetary collisions. For both HR8799b and 
2M1207, iBarman et al.l (l2011al Jbl) point out that the thick clouds (reminiscent of L-dwarfs) 
extending across the photosphere, effectively accounting for the observed red colors, and the 
non-equilibrium chemistry consistent with the observed CO/CH4 ratio are equally important 
but only in the presence of low surface gravity. 



Further observations of the HR8799 planets (jSkemer et al.ll2012l ) also indicate that all 
four planets are brighter than expected at 3.3 fim compared to equilibrium models that pos- 
tulate significant methane absorption and dimming at this wavelength. The authors report 
that the observations are inconsistent with models with decreased CH4, thick clouds and 
non-equilibrium che mistry. Their 3 .3 nm photometry of HR8799b is not co nsistent with 



the be st-fit model of IBarman et al.l (1201 laf) and is twi ce the value obtained by ICurrie et al. 
(I2OIII). Using the rn odels of lMadhusudhan et al.l ( I2OIII ) and adopting a "patchy" cloud cover, 
Skemer et al.l (120121 ) find a best-fit model (to all their p hotometry data excep t M-b and) with 
a two-component atmosphere of 93% Teff = 700 K (IMadhusudhan et al.l 1201 ll )('A'-type 
clouds) and 7% Teff = 1400 K ('AE'-type clouds) for HR8799b. Similar hybrid model at- 
mospheres provide better fit to HR8799 c, d, and e compared to thick clouds/non-equilibrium 
chemistry models that again fail to reproduce the 3.3 /xm colors. The four H R8799 planets 
have d ifferent effective temperatures but, intriguingly, similar colors which led lSkemer et al. 
(I2OI2I ) to suggest that their atmospheres should have similar properties and may indeed be 
an evidence for their patchy appearance. 



Low-level photometric variability due to rotation of 



by multiple groups in the IRAC 4.5 fi m and 8 fim bands ( Morales- Calderon et al.ll2006f). and 



at sh orter near- infrared wavelengths (IClarke et al.l 12008 



and T dwarfs have been rep orted 



Artigau et al. 



2009 



Radigan et al. 



2OI2I ). Ongoing large Spitzer and Hubble Space Telescope programs are identifying a larger 
number of previously unknown varying brown dwarfs. 



5.2. Photometric and Spectroscopic Variability of Directly-Imaged Giant 

Planets 

We showed that rotation periods can be measured with existing adaptive optics systems 
(Section [3^ if the modulation is as prominent as our models predict (~ 20%). The detection 
of variability with lower amplitudes will be possible with next-generation adaptive optics 
systems. One such system (LBT/AO), for example, already provides such a capability 



( ISkemer et al.l |2 0121 ). Rotational variability can provide three different types of information 
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for directly imaged exoplanets: 1) single-band photometry can determine the rotational 
period without the ambiguity of inclination; 2) multi-band photometry can constrain the 
heterogeneity of the surface and the relative colors of the features, and 3) spectral mapping 
can provide spectrally and spatially resolved maps of the ultra cool atmospheres, offering 
detailed insights into the atmospheric structures of giant exoplanets. In the following we 
will briefly discuss these different measurements and their observational requirements as 
predicted by our simple model. 

We showed that rotation periods can already be measured with existing adaptive optics 
systems (Section l3.2p . We find that for planets broadly similar to those in the HR 8799 system 
the amplitude of predicted variations peaks in the J-band, yet the optimal wavelength for 
observations is the Ks-band, where AO performance is superior. 

Next-generation AO systems on large telescopes will provide greatly improved AO per- 
formance at shorter wavelengths. Our models predict that these systems will be able to 
detect rotational variations at multiple wavelengths, determining cloud covering fraction 
and providing color information on the surface features (see Section 1X5]) . 

We also briefly explored spectral mapping of exoplanets on the example of JWST. These 
observations will provide spatially and spectrally resolved maps over a broad wavelength 
range and unaffected by telluric absorption. Our models predict that such observations 
both with JWST and with instruments on 30m telescopes will be limited by instrumental 
stability and not by contrast. This fact highlights the importance for sub-percent-level flux 
calibration techniques. 

We also use our models to evaluate the factors that determine the detectability of 
rotational variations. We find that temperature, spatial and size distributions of the spots 
and the effective temperature of the planet play the most important roles. For the spot 
distributions we modeled (Fig. [1]), changing the inclination to +30° does not affect the 
amplitude of the lightcurves compared to 0° inclination but slightly changes their shape. On 
the contrary, an inclination of —50° emphasizes the contribution from the two cold spots 
in the right panel of Figure [T] and changes not only the shape of the lightcurves but also 
the maximum amplitudes (see Figured!]) - they decrease by 12%, 9.5% and 10% in the J-, 
H- and Ks-bands respectively compared to the maximum amplitudes for the 0° inclination 
scenario seen in Figure |H The importance of the inclination of the planet also depends on 
the spatial distribution of the spots - if they are distributed mostly around the equatorial 
region, highly-inclined planets will show lower level variability than planets seen from their 
equatorial plane. 

We note the following major difference between the Clear and Cloudy B06 models we 
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use that have major effects on the results presented in this paper: a) The ffux densities 
of the Cloudy (yellow, Figure [3]) and Clear (green, Figure |3]) B06 models of 1,000 K are 
similar in both H and K-bands but differ significantly in the J-band, where the Clear model 



spectrum i s notably higher. T his is the well-known J-brightening effect (iBurgasser et al. 



2000l . I2OO2I : iLeggett et al.ll2000l ) caused by the appearance of cloud- free regions at the L/T 
transition; b) The J-band is stronger than the H-band for the 1,000 K Cloudy (yellow. 
Figure [3]) but the two bands have similar strength for the 1,200 K Cloudy (light orange. 
Figure E]) models. The Ks-band is significantly weaker for both temperatures; and c) the 
1,400 K Cloudy (dark orange. Figure |3]) model spectrum peaks in the H-band and the flux 
densities in the J and Ks-bands are somewhat similar. Therefore, we can expect that a 
heterogenous surface consisting of multiple patches (spots, holes and/or clouds) described 
by different spectra will have significantly different observational signatures compared to a 
homogeneous, one-surface-type atmosphere. 

The distribution and sizes of cloud structures on directly imaged giant exoplanets will 
not be known a priori but constraints can be obtained by the model presented here. In Section 
3 we showed that even a moderate spot total covering factor of 10% and spot distribution 
similar to that of Jupiter with a dominant giant spot can produce up to 20% photometric 
variability in J, H and Ks bands for a Cloudy surface with Tg// of 800 K, 1,000 K and 
1,200 K, as seen from Figure [71 Also, a giant spot with a temperature difference as small 
as 6T = 100 K compared to the effective temperature can cause photometric modulations 
of up to 20% in Ks-band for an effective temperature of T^jf = 800 K. We also argued that 
the most appropriate filter, where the largest simulated photometric variability consistently 
occurs depends on the different surface types present on the planet. As discussed in Section 
13. H for our Model B2 and for the map in Figure[T]the filter with the largest variations is Ks for 
Teff lower than 1,000 K and J for Teff ranging from 1,000 K to 1,400 K. However, for Model 
A2 the largest photometric variations occur in Ks for all Teff from TOOK to 1,400 K. This 
suggests that simply detecting the wavelength at which the largest photometric variations 
occur will already be a strong indicator of the relative contributions of cloudy and cloud-free 
regions to the planet's atmosphere and also of the temperatures of these regions. 



Barman et al.l (l2011a] ) report a weak CH4 absorption in both H- and K-bands and a 
triangular shape in the H-band spectrum for the case of HR8799b, indicating low surface 
gravity which may promote more efficient vertical mixing, deviation from chemical equilib- 
rium and indicate a young age. As noted by the authors, such low surface gravity would also 
imply that the condensation curve crosses the T-P curve near the photosphere and clouds 
can form in the deeper, photospheric depths, suggesting that the giant young planets can be 
cloudy and cool. For the same effective temperature and surface gravity, the thickness of the 
cloud layer is significantly smaller compared to the case of chemical equilibrium but still suf- 
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ficient enough to produce the observed colors, much redder than those of a c loud-free brown 



dwar f. Cloudy, cool, low-gravity brown dwarfs have indeed been reported (jStephens et al. 



20091 ) at the L/T transition, pushing the transition temperature from 1,300 K (at log(g) of 
5) down to 1,100 K (at log(g) of 4.5), supporting e arlier evidence for such g ravity-transition 
regime connection (IMetchev fc Hillenbrandl l2006t ISaumon fc Marleyll2008l ). The rotational 
phase mapping proposed here can provide critical insights into the structure and distribu- 
tion of cloud layers and can test cloud scale height and compositional variations in directly 
imaged giant planets (Section [3]). 

Our model predicts that the detection of photometric and spectroscopic variations with 
current instrumentation are limited by a combination of contrast and stability (section |3]), 
but next-generation extreme AO systems will provide superior contrasts that will no longer 
be the limiting factor. The sensitivity limits of the current state-of-the-art instruments like 
VLT/NACO limit the capabilities of these instruments to measuring rotation periods. With 
next-generation instruments like GPI, SPHERE, or extreme AO systems on GMT, TMT 
or E-ELT will be able to explore the details of the cloud asymmetry through the lower- 
amplitude (down to an assumed photometric precision of 1%) rotational modulations and 
study the J— Ks and H— Ks colors of spots covering the same fraction (~ 1%) of the surface 
area while JWST may push this limit down to the level of ~ 10~^. The real power of the 
future very large-aperture telescopes, both ground- and space-based like GMT, TMT, E- 
ELT, JWST and possibly ATLAST, will be in the spectral modulation domain, where they 
will open up the possibility to study in detail the cloud colors, composition and/or spectra, 
weather patterns and even rings and satellites. 

T he Great Red Spot on J upiter has already been shown to be a very dynamic fea- 



ture - lAsay-Davis et al.l (120091 ) showed that the spot has shrunk by 15% over a period of 
10 years. Such behavior needs to be taken into account when studying the atmospheres of 
directly-imaged giant exoplanets - variations on timescales different than the rotation period 
can indicate evolving atmospheric patterns or even storms. If the atmospheric patterns of 
these planets do indeed change on timescales shorter than the rotation period, the periodic 
photometric a nd/or spectrosc opic modulations will be significantly modified or even com- 
pletely erased ( lGoldmanll2005l ) . Detections of photometric and/or spectroscopic evolution on 
timescales different from the rotation period will be a further step toward understanding the 
atmosphere of giant exoplanets. Such measurements, possibly achievable by JWST as dis- 
cussed in Section 13. 5[ w ill require very high- c adence observation s and could be an indication 
of differential rotation (lArtigau et al.l l2009l : iMarley et al.l 120101 ). We have not commented 
on the possible presence of symmetric belts or bands in the atmosphere of these planets as 
they will not cause photometric or spectroscopic variability and cannot be constrained by 
the model presented here. 
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6. Conclusions 

We used simple models to explore light curve inversion as a tool to probe the surface 
brightness distribution and cloud properties of directly imaged exoplanets. We constructed 
models using state-of-the-art spectral libraries and simulated observations with current and 
future instruments. The key results of this study are as follows: 

i) Current AO systems on large telescopes are capable of detecting large photometric 
variations and thus measure rotation periods. Cadence is limited by the slowly varying PSF 
systematics. 

ii) Next-generation extreme AO systems will be capable of detecting even low-level 
(~1%) variations, enabling the study of cloud heterogeneity and its wavelength dependence. 
These systems will also enable the study of the cloud cover evolution. 

iii) Extremely large telescopes and JWST will provide spectral mapping data for clouds, 
allowing detailed composition maps of the cloud cover and the abundance of gas-phase 
absorbers. These setups will also provide a much higher cadence. 

iv) For objects at the L/T transition, a B2 model with T^ff = lOOOi^T and a 1200ir hot 
spot predicts the largest photometric variations to be in the J-band; for cooler sources or 
sources with very cool surface features the ideal wavelength gradually shifts toward 3 /xm. 

v) We demonstrated that simulated data can be inverted to a correct, low-resolution 
one-dimensional map of a giant planet. 

The observations proposed here will allow detailed studies of the structure and com- 
position of condensate cloud covers in directly imaged exoplanets and provide otherwise 
inaccessible insights into the atmospheric circulation of these exciting objects. 

7. Acknowledgments 

A Director's Discretionary research grant by the Space Telescope Science Institute was 
essential for starting this program. We are grateful to Markus Kasper for valuable discus- 
sions. We would also like like to acknowledge Justin Rogers for encouraging discussions and 
for his careful review of this paper. 



-27- 



REFERENCES 

AUard, F., Hauschildt, P. H., Alexander, D. R., Tamanai, A., & Schweitzer, A. 2001, ApJ, 
556, 357 

AUard, F., Homeier, D., & Frcytag, B. 2011, in Astronomical Society of the Pacific Confer- 
ence Series, Vol. 448, 16th Cambridge Workshop on Cool Stars, Stellar Systems, and 
the Sun, ed. C. Johns-KruU, M. K. Browning, & A. A. West, 91 

Apai, D., Janson, M., Moro-Martm, A., Meyer, M. R., Mamajek, E. E., Masciadri, E., 
Henning, T., Pascucci, I., Kim, J. S., Hillenbrand, L. A., Kasper, M., & Biller, B. 
2008, ApJ, 672, 1196 

Apai, D. et al. 2012, in prep. 

Artigau, E., Bouchard, S., Doyon, R., & Lafreniere, D. 2009, ApJ, 701, 1534 

Asay-Davis, X. S., Marcus, P. S., Wong, M. H., & de Pater, I. 2009, Icarus, 203, 164 

Bailer- Jones, C. A. L. 2008, MNRAS, 384, 1145 

Bailer- Jones, C. A. L. & Mundt, R. 1999, A&A, 348, 800 

— . 2001, A&A, 367, 218 

Baines, K. H., Drossart, P., Momary, T. W., Formisano, V., Griffith, C, Bellucci, G., 
Bibring, J. P., Brown, R. H., Buratti, B. J., Capaccioni, F., Cerroni, P., Clark, R. N., 
Coradini, A., Combes, M., Cruikshank, D. P., Jaumann, R., Langevin, Y., Matson, 
D. L., McCord, T. B., Mennella, V., Nelson, R. M., Nicholson, P. D., Sicardy, B., & 
Sotin, C. 2005, Earth Moon and Planets, 96, 119 

Baraffe, 1., Chabrier, G., Barman, T. S., AUard, F., & Hauschildt, P. H. 2003, A&A, 402, 
701 

Barman, T. S., Macintosh, B., Konopacky, Q. M., & Marois, C. 2011a, ApJ, 733, 65 
— . 2011b, ApJ, 735, L39 

Barrado-Izagirre, N., Perez-Hoyos, S., & Sanchez-Lavega, A. 2009, Icarus, 202, 181 

Bonnefoy, M., Lagrange, A.-M., Boccalctti, A., Chauvin, G., Apai. D., Allard, F., Ehrenre- 
ich, D., Girard, J. H. V., Mouillet, D., Rouan, D., Gratadour, D., & Kasper, M. 2011, 
A&A, 528, L15 



-28- 



Burgasser, A. J. 2009, in American Institute of Physics Conference Series, Vol. 1094, 15th 
Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, ed. E. Stempels, 
501-504 

Burgasser, A. J., Marley, M. S., Ackerman, A. S., Saumon, D., Lodders, K., Dahn, C. C, 
Harris, H. C, & Kirkpatrick, J. D. 2002, ApJ, 571, L151 

Burgasser, A. J., Wilson, J. C, Kirkpatrick, J. D., Skrutskic, M. F., Colonno, M. R., Enos, 
A. T., Smith, J. D., Henderson, C. P., Gizis, J. E., Brown, M. E., & Houck, J. R. 
2000, AJ, 120, 1100 

Burrows, A. 2009, in Astronomical Society of the Pacific Conference Series, Vol. 414, Cosmic 
Dust - Near and Far, ed. T. Henning, E. Griin, & J. Steinacker, 115 

Burrows, A., Sudarsky, D., & Hubeny, I. 2006, ApJ, 640, 1063 

Chabrier, C, Baraffe, 1., AUard, F., & Hauschildt, P. 2000, ApJ, 542, 464 

Chauvin, C, Lagrange, A.-M., Dumas, C, Zuckerman, B., Mouillet, D., Song, I., Beuzit, 
J.-L., & Lowrance, P. 2004, A&A, 425, L29 

Clampin, M. 2008, Advances in Space Research, 41, 1983 

Clarke, F. J., Hodgkin, S. T., Oppenheimer, B. R., Robertson, J., & Haubois, X. 2008, 
MNRAS, 386, 2009 

Cowan, N. B., Agol, E., Meadows, V. S., Robinson, T., Livengood, T. A., Deming, D., Lisse, 
C. M., A'Hearn, M. F., Wellnitz, D. D., Seager, S., Charbonneau, D., & the EPOXI 
Team. 2009, ApJ, 700, 915 

Currie, T., Burrows, A., Itoh, Y., Matsumura, S., Fukagawa, M., Apai, D., Madhusudhan, 
N., Hinz, P. M., Rodigas, T. J., Kasper, M., Pyo, T.-S., & Ogino, S. 2011, ApJ, 729, 
128 

Cushing, M. C, RoeUig, T. L., Marley, M. S., Saumon, D., Leggett, S. K., Kirkpatrick, 
J. D., Wilson, J. C, Sloan, G. C, Mainzer, A. K., Van Cleve, J. E., & Houck, J. R. 
2006, ApJ, 648, 614 

Dekany, R. G., Britton, M. C, Gavel, D. T., EUerbroek, B. L., Herriot, G., Max, C. E., 
& Veran, J.-P. 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 5490, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, ed. D. Bonaccini Calia, B. L. EUerbroek, & R. Ragazzoni, 
879-890 



-29- 



Desert, J.-M., Bean, J., Miller-Ricci Kempton, E., Berta, Z. K., Charbonneau, D., Irwin, J., 
Fortney, J., Burke, C. J., & Nutzman, P. 2011, ApJ, 731, L40 

Dyudina, U. A., IngersoU, A. P., Danielson, G. E., Baines, K. H., Carlson, R. W., The 
Galileo Nims, & SSI Teams. 2001, Icarus, 150, 219 

Fortney, J. J., Marley, M. S., Saumon, D., & Lodders, K. 2008, ApJ, 683, 1104 

Gardner, J. P., Mather, J. C., Clampin, M., Doyon, R., Greenhouse, M. A., Hammel, H. B., 
Hutchings, J. B., Jakobsen, P., Lilly, S. J., Long, K. S., Lunine, J. I., McCaughrean, 
M. J., Mountain, M., Nella, J., Rieke, G. H., Rieke, M. J., Rix, H.-W., Smith, E. P., 
Sonneborn, G., Stiavelh, M., Stockman, H. S., Windhorst, R. A., & Wright, G. S. 
2006, Space Sci. Rev., 123, 485 

Gelino, C. & Marley, M. 2000, in Astronomical Society of the Pacific Conference Series, Vol. 
212, From Giant Planets to Cool Stars, ed. C. A. Griffith & M. S. Marley, 322 

GMT. 2006, GMT Science Case, http://www.gmto.org/sciencecase.html 

Goldman, B. 2005, Astronomische Nachrichten, 326, 1059 

Green, J. J., Beichman, C, Basinger, S. A., Horner, S., Meyer, M., Redding, D. C, Rieke, M., 
& Trauger, J. T. 2005, in Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 5905, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, ed. D. R. Coulter, 185-195 

Herriot, G., Hickson, P., EUerboek, B. L., Andersen, D. A., Davidge, T., Erickson, D. A., 
Powell, I. P., Clare, R., Gilles, L., Boyer, C, Smith, M., Saddlemyer, L., & Veran, 
J.-P. 2006, in Society of Photo-Optical Instrumentation Engineers (SPIE) Confer- 
ence Series, Vol. 6272, Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series 

Hinz, P. M., Rodigas, T. J., Kenworthy, M. A., Sivanandam, S., Heinze, A. N., Mamajek, 
E. E., & Meyer, M. R. 2010, ApJ, 716, 417 

Irwin, P. G. J. 2003, Giant planets of our solar system : atmospheres compositions, and 
structure 

Jaffe, D. T., Dcpoy, D. L., Fabricant, D. G., Hinz, P. M., Jacoby, G., Johns, M., McCarthy, 
P., McGregor, P. J., Shectman, S., & Szentgyorgyi, A. 2010, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7735, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 



-30- 



Johns, M. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 6986, Society of Photo-Optical Instrumentation Engineers (SPIE) Con- 
ference Series 

Kasper, M., Amico, P., Pompei, E., Ageorges, N., Apai, D., Argomedo, J., Kornweibel, N., 
& Lidman, C. 2009, The Messenger, 137, 8 

Kasper, M., Apai, D., Janson, M., & Brandner, W. 2007, A&A, 472, 321 

Kasper, M., Beuzit, J.-L., Verinaud, C, Gratton, R. G., Kerber, P., Yaitskova, N., Boc- 
caletti, A., Thatte, N., Schmid, H. M., Keller, C, Baudoz, P., Abe, L., Allcr- 
Carpentier, E., Antichi, J., Bonavita, M., Dohlen, K., Fedrigo, E., Hanenburg, H., 
Hubin, N., Jagcr, R., Korkiakoski, V., Martinez, P., Mesa, D., Preis, O., Rabou, 
P., Roclfscma, R., Salter, G., Tccza, M., & Venema, L. 2010, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7735, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 

Kasper, M. E., Beuzit, J.-L., Verinaud, C, Yaitskova, N., Baudoz, P., Boccaletti, A., Grat- 
ton, R. G., Hubin, N., Kerber, F., Roelfsema, R., Schmid, H. M., Thatte, N. A., 
Dohlen, K., Feldt, M., Venema, L., & Wolf, S. 2008, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 7015, Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series 

Knapp, G. R., Leggett, S. K., Fan, X., Marley, M. S., Geballe, T. R., Golimowski, D. A., 
Finkbeiner, D., Gunn, J. E., Hennawi, J., Ivezic, Z., Lupton, R. H., Schlegel, D. J., 
Strauss, M. A., Tsvetanov, Z. I., Chiu, K., Hoversten, E. A., Glazebrook, K., Zheng, 
W., Hendrickson, M., Williams, C. C, Uomoto, A., Vrba, F. J., Henden, A. A., 
Luginbuhl, C. B., Guetter, H. H., Munn, J. A., Canzian, B., Schneider, D. P., & 
Brinkmann, J. 2004, AJ, 127, 3553 

Lafreniere, D., Marois, C, Doyon, R., Nadeau, D., & Artigau, E. 2007, ApJ, 660, 770 

Lagage, P. O. & European MIRI Team. 2010, in In the Spirit of Lyot 2010 

Lagrange, A.-M., Bonncfoy, M., Chauvin, G., Apai, D., Ehrenreich, D., Boccaletti, A., 
Gratadour, D., Rouan, D., Mouillet, D., Lacour, S., & Kasper, M. 2010, Science, 329, 
57 

Lagrange, A.-M., Gratadour, D., Chauvin, G., Fusco, T., Ehrenreich, D., Mouillet, D., 
Rousset, G., Rouan, D., AUard, F., Gendron, E., Charton, J., Mugnier, L., Rabou, 
P., Montri, J., & Lacombe, F. 2009, A&A, 493, L21 



-31 - 



Leggett, S. K., Geballe, T. R., Fan, X., Schneider, D. P., Gunn, J. E., Lupton, R. H., 
Knapp, G. R., Strauss, M. A., McDaniel, A., Golimowski, D. A., Henry, T. J., Peng, 
E., Tsvetanov, Z. I., Uomoto, A., Zheng, W., HiU, G. J., Ramsey, L. W., Anderson, 
S. P., Annis, J. A., Bahcall, N. A., Brinkmann, J., Chen, B., Csabai, I., Pukugita, M., 
Hennessy, G. S., Hindsley, R. B., Ivezic, Z., Lamb, D. Q., Munn, J. A., Pier, J. R., 
Schlegel, D. J., Smith, J. A., Stoughton, C., Thakar, A. R., & York, D. G. 2000, ApJ, 
536, L35 

Liske, J. et al 2011, The E-ELT Design Reference Mission, E-TRE-ESO-080-0717, Issue 2 

Macintosh, B., Troy, M., Doyon, R., Graham, J., Baker, K., Bauman, B., Marois, C., Palmer, 
D., PhiUion, D., Poyneer, L., Crossfield, I., Dumont, P., Levine, B. M., Shao, M., Ser- 
abyn, G., Shelton, C., Vasisht, G., Wallace, J. K., Lavigne, J. -P., Valee, P., Rowlands, 

N., Tam, K., & Hackett, D. 2006, in Society of Photo-Optical Instrumentation Engi- 
neers (SPIE) Conference Series, Vol. 6272, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series 

Macintosh, B. A., Graham, J. R., Palmer, D. W., Doyon, R., Dunn, J., Gavel, D. T., 
Larkin, J., Oppenheimer, B.. Saddlemyer, L., Sivaramakrishnan, A., Wallace, J. K., 
Bauman, B., Erickson, D. A., Marois, C, Poyneer, L. A., & Soummer, R. 2008, in 
Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 
7015, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series 

Madhusudhan, N., Burrows, A., & Currie, T. 2011, ApJ, 737, 34 

Marley, M. S., Saumon, D., & Goldblatt, C. 2010, ApJ, 723, L117 

Marois, C., Lafreniere, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006a, ApJ, 641, 556 

Marois, C, Lafreniere, D., Macintosh, B., & Doyon, R. 2006b, ApJ, 647, 612 

Marois, C, Macintosh, B., Barman, T., Zuckerman, B., Song, I., Patience, J., Lafreniere, 
D., & Doyon, R. 2008, Science, 322, 1348 

Marois, C, Zuckerman, B., Konopacky, Q. M., Macintosh, B., & Barman, T. 2010, Nature, 
468, 1080 

McCarthy, Jr., D. W., Ge, J., Hinz, J. L., Finn, R. A., & de Jong, R. S. 2001, PASP, 113, 
353 

McLean, 1. S., ed. 1997, Electronic imaging in astronomy. Detectors and instrumentation 



-32 - 



Mesa, D., Gratton, R., Berton, A., Antichi, J., Verinaud, C, Boccaletti, A., Kasper, M., 
Claudi, R. U., Desidera, S., Giro, E., Beuzit, J.-L., Dohlen, K., Feldt, M., Mouillet, 
D., Ghauvin, G., & Vigan, A. 2011, A&A, 529, A131 

Metchev, S. A. & Hillenbrand, L. A. 2006, ApJ, 651, 1166 

Miller-Ricci Kempton, E., Zahnle, K., & Fortney, J. J. 2011, ArXiv e-prints 

Mohanty, S., Jayawardhana, R., Huelamo, N., & Mamajek, E. 2007, ApJ, 657, 1064 

Morales-Caldcron, M., Stauffcr, J. R., Kirkpatrick, J. D., Carey, S., Gclino, C. R., Barrado 
y Navascues, D., RebuU, L., Lowrance, P., Marley, M. S., Charbonneau, D., Patten, 
B. M., Megeath, S. T., & Buzasi, D. 2006, ApJ, 653, 1454 

Oppenheimer, B. R., Digby, A. P., Ncwburgh, L., Brenner, D., Shara, M., Mey, J., Mandev- 
ille, C., Makidon, R. B., Sivaramakrishnan, A., Soummer, R., Graham, J. R., Kalas, 
P., Pcrrin, M. D., Roberts, Jr., L. C., Kuhn, J. R., Whitman, K., & Lloyd, J. P. 2004, 
in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 
5490, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 
ed. D. Bonaccini Calia, B. L. EUerbroek, & R. Ragazzoni, 433-442 

Orton, G., Ortiz, J. L., Baincs, K., Bjorakcr, G., Carsenty, U., Colas, P., Dayal, A., Deming, 
D., Drossart, P., Frappa, E., Friedson, J., Goguen, J., Golisch, W., Griep, D., Her- 
nandez, C., Hoffmann, W., Jennings, D., Kaminski, C, Kuhn, J., Laques, P., Limaye, 
S., Lin, H., Lecacheux, J., Martin, T., McCabe, G., Momary, T., Parker, D., Puetter, 
R., Ressler, M., Reyes, G., Sada, P., Spencer, J., Spitale, J., Stewart, S., Varsik, J., 
Warell, J., Wild, W., Yanamandr a- Fisher, P., Fazio, G., Hora, J., & Deutsch, L. 1996, 
Science, 272, 839 

Patience, J., King, R. R., de Rosa, R. J., & Marois, C. 2010, A&A, 517, A76 

Postman, M., Brown, T., Sembach, K., Giavalisco, M., Traub, W., Stapclfeldt, K., Calzetti, 
D., Oegerle, W., Rich, R. M., Stahl, H. P., Tumlinson, J., Mountain, M., Soummer, 
R., & Hyde, T. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 7731, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series 

Quanz, S. P., Meyer, M. R., Kenworthy, M. A., Girard, J. H. V., Kasper, M., Lagrange, 
A.-M., Apai, D., Boccaletti, A., Bonnefoy, M., Ghauvin, G., Hinz, P. M., & Lenzen, 
R. 2010, ApJ, 722, L49 



-33- 



Radigan, J., Jayawardhana, R., Lafreniere, D., Artigau, E., Marley, M., & Saumon, D. 2012, 
ApJ, 750, 105 

Rogers, J. H. 2008, Journal of the British Astronomical Association, 118, 14 

Rogers, L. A. & Seager, S. 2010, ApJ, 716, 1208 

Russell, H. N. 1906, ApJ, 24, 1 

Saumon, D. & Marley, M. S. 2008, ApJ, 689, 1327 

Seager, S. & Deming, D. 2010, ARA&A, 48, 631 

SeifF, A., Kirk, D. B., Knight, T. C. D., Young, R. E., Mihalov, J. D., Young, L. A., Milos, 
F. S., Schubert, G., Blanchard, R. C, & Atkinson, D. 1998, J. Geophys. Res., 1032, 
22857 

Showman, A. P., Cho, J. Y.-K., & Menou, K. Atmospheric Circulation of Exoplanets, ed. 
S. Seager, 471-516 

Skemer, A. J., Close, L. M., Szucs, L., Apai, D., Pascucci, I., & Biller, B. A. 2011, ApJ, 732, 
107 

Skemer, A. J., Hinz, P. M., Esposito, S., Burrows, A., Leisenring, J., Skrutskie, M., Desidera, 
S., Mesa, D., Arcidiacono, C, Mannucci, F., Rodigas, T. J., Close, L., McCarthy, D., 
Kulesa, C, Agapito, C, Apai, D., Argomedo, J., Bailey, V., Boutsia, K., Briguglio, 
R., Brusa, C, Busoni, L., Claudi, R., Eisner, J., Fini, L., FoUette, K. B., Garnavich, 
P., Gratton, R., Guerra, J. C, Hill, J. M., Hoffmann, W. F., Jones, T., Krejny, M., 
Males, J., Masciadri, E., Meyer, M. R., Miller, D. L., Morzinski, K., Nelson, M., 
Pinna, E., Puglisi, A., Quanz, S. P., Quiros-Pacheco, F., Riccardi, A., Stefanini, P., 
Vaitheeswaran, V., Wilson, J. C, & Xompero, M. 2012, ArXiv e-prints 

Snyder, J. P. Map Projections - A Working Manual, ed. Washington, D.C.: US Government 
Printing Office, 145-153 

Stephens, D. C, Leggett, S. K., Gushing, M. C., Marley, M. S., Saumon, D., Geballe, T. R., 
Golimowski, D. A., Fan, X., & Noll, K. S. 2009, ApJ, 702, 154 

Stiavelh, M. et al. 2008, "JWST Primer", Version 2.0, (Baltimore: STScI) 

Tsuji, T. & Nakajima, T. 2003, ApJ, 585, L151 

Vasavada, A. R. & Showman, A. P. 2005, Reports on Progress in Physics, 68, 1935 



-34- 



Verinaud, C, Kasper, M., Beuzit, J.-L., Gratton, R. G., Mesa, D., AUer-Carpentier, E., 

Fedrigo, E., Abe, L., Baudoz, P., Boccaletti, A., Bonavita, M., Dohlcn, K., Hubin, 
N., Kerber, P., Korkiakoski, V., Antichi, J., Martinez, P., Rabou, P., Roelfsema, R., 
Schmid, H. M., Thatte, N., Salter, G., Tecza, M., Venema, L., Hanenburg, H., Jager, 
R., Yaitskova, N., Preis, O., Orecchia, M., & Stadler, E. 2010, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7736, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 

Vigan, A., Moutou, C, Langlois, M., Allard, F., Boccaletti, A., Carbillet, M., Mouillet, D., 
& Smith, 1. 2010, MNRAS, 407, 71 

Westphal, J. A., Matthews, K., & Terrile, R. J. 1974, ApJ, 188, Llll 



This preprint was prepared with the A AS macros v5.2. 



- 35 - 



Orthographic projection of the model mop 




Front Bock 
Orange = 1000K, Yellow = 1 200K, Red = 800K 



Fig. 1. — Projected model surface map of a giant planet with randomly distributed spots ac- 
cording to a power law with an index of -2.0 (such that the largest five clouds are responsible 
for 85% of the signal). Different colors correspond to different surface types temperatures, 
as indicated in the figure. This model has a 10% total spot-covering factor (the giant spot 
covers 5% of the total surface area), 3 different temperatures, and inclination of 0°. The left 
disk represents the front side of the planet facing the observer, the right disk - the back. 
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Fig. 2. — The visible surface fraction covered by the 5 largest spots as a function of rotational 
phase for the map shown in Figure [TJ 
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Model spectra from Burrows and AMES — Cond 
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Fig. 3. — Sample spectra for temperatures of 800K, lOOOK and 1200K and log(g)=4.5 from 
the Clear Burrows (B06, top) and AMES-Cond (AOl, lower) models. 
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Model lightcurves for JHKL filters 
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Fig. 4. — Normalized model lightcurves as a function of rotation phase for the map in 
Figure [H using Model B2 (clear spots on a cloudy surface, with different temperatures) and 
assuming Tg// = lOOO-ft' and the Giant Spot at 1200K. As discussed in the text, a suite 
of next-generation instruments, dedicated to direct imaging of giant exoplanets, can detect 
such modulations and measure the rotation periods and/or cloud-coverage of the planet. 
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Spectral Voriations from o Rotating Giant Planet with 3 surface types 
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Fig. 5. — Predicted spectral modulations (normalized) for JWST/NIRSPEC and 
JWST/MIRI (with a resolution of 100) as a function of the rotation phase for the map 
in Figure [H using a model with clear spots on a cloudy surface (B2 from Table [3]) for a 
planet with T^fj =1,000 K and a giant hot spot of 1,200 K. The top panel shows the visible 
hemisphere as a function of orbital phase (shown at the bottom). Absorption bands for the 
respective molecules are marked on the right. 
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Rototion Phose 

Fig. 6. — Same as Figure Obut for different properties of the spots (temperature and surface 
type). Top panel is for Model Al, middle panel is for Model A2 and lower panel is for Model 
C. All models are for a T^ff = lOOOif . 
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Amplitude of lightcurve variations 
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Fig. 7. — Normalized amplitudes of lightcurve variations as a function of the temperature of 
the giant spot for the map in Figured] The top panel is for T^fj = 800K, the middle panel 
is for Teff = lOOOif and the lower panel - for T^ff = 12007^'. All three panels are for Model 
B2 (Table [3]). The largest variations occur in different filters depending on the temperature 
of the giant spot, suggesting that the observations should be carefully tailored to the specific 
target. 
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Simulated Observations 
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Fig. 8. — Simulated observations of a Tg// = lOOO/T planet with a giant hot spot of 1200K 
(Model B2) for an 8-m class telescope with AO (upper panel) and with Extreme AO (mid- 
dle panel), and for a 30-m class telescope with Extreme AO (lower panel). The red and 
yellow line are the theoretical lightcurves from Figure HI the black diamonds are simulated 
photometry; the error bars represent a photometric precision of 10% in the upper panel and 
1% in the middle and lower panels. We assume a rotation period of 4 hours. The temporal 
sampling is 15 min cadence for the 8-m AO, 8 min for the 8-m ExAO and 10 min cadence 
for the 30-m aperture. 
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SImuloted Observotlons from JWST 
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Fig. 9. — Same as the previous figure, but for simulated observations from JWST/NIRCAM 
F115W witli a cadence of 15 min and photometric precision of ~ 10~^. 
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Principal Eigenvalues by Weight 
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Fig. 10.— Principal eigenvalues by weight: 0.998, 2.2xl0-^ 2.2 xlO'^ 6.3x10"^ and 
3.1xl0~*. The two largest values correspond to the two eigencolors representing the two 
different spot temperatures we assumed in our model. The results shows that the simulated 
photometric variations are indeed caused by only two colors (different from the mean color, 
which is a manifest of the effective temperature of the planet), or precisely the number of 
"extra colors" in our input map. The other three eigenvalues are zero to the precision we 
used. We used our own PCA routine, but the built-in IDL PCA routine also gives zeros for 
these three. 
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Fig. 11. — Phase variations of the primary (top panel) and secondary (lower panel) eigen- 
colors as inverted by PCA from the simulated photometry from Figure HI Rotation phase 
of zero is defined as the left panel on [TJ A maximum in the relative color outlines a rota- 
tion phase at which the most amount (largest area covered) of that respective eigencolor is 
present on the side of the planet facing the observer and a minimum - the least. The primary 
eigencolor practically mimics the input (orange) J-band lightcurve from Figure IH while the 
secondary has a completely different behavior, hard to notice by eye from the lightcurve. 
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Fig. 12. — MoUweide projection of the input map (upper panel: orange = lOOOi^, yellow = 
1200i^', black = 8007^) and of the PCA-inverted, longitudinal distribution of the primary 
eigencolor (lower panel: normalized percentage contribution). The inversion method used 
successfully recovers the longitudinal position (there is no latitudinal resolution) of the giant 
hot spot at +135° and, to a somewhat lesser degree, that of the two smaller spots at ~ —110°. 
The secondary eigencolor is masked out (as black) in the upper panel as it is orthogonal to 
the primary and does not contribute to the variations of the primary. 
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Longitudinol Map of Secondory Eigencolor (Red) 
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Fig. 13. — The same as the previous figure, but for the secondary eigencolor (upper panel: 
orange = lOOOi^, black = 12007^, red = SOOi^'). As in the case for the primary color, the 
inversion successfully recovers the longitudinal position of the two groups of cold spots, at 
~ —60° and ~ +60° respectively. The primary eigencolor is masked out (as black) in the 
upper panel. 
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Orthogrophic projection, inclination -50 degrees 
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Fig. 14. — Top panel: the same as Figured] but for an inclination of —50°. Lower panel: 
the same as Figure H] but again for an inclination of —50°. The maximum amplitudes in all 
filters are significantly smaller (by as much as 12% in J-band) compared to an inclination of 
0°. The shapes of the lightcurves change as well (compared toH]), emphasizing the lessened 
contribution from the giant hot spots. 
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Telescope 


Instrument 


Wavelength [/im] 


Reference 


8-m class ExAO 


SPHERE/GPI/PISCES 


1 to 2.5 


1 


30-m class ExAO 


GMT/TMT/ELT 


1 to 5 


2 


JWST 


NIRCAM/MIRI 


1 to 27 


3 


ATLAST 




0.1 to 2.4 


4 



Table 1: List o: 
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Postman et al. 
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20101 ): 
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Ta ble 2: Sensitiv i ty lim i ts for the instruments we exp l ore at the respectiv e an gular separation: 



Kasper et all (12007 . 



(120081 



20091) 



McCarthv et al.l (120011): 



Kasper et all J2008I . boioh [7] 



[2] IVigan et al.l (120101) : 



Skemer et al. 



Green et al 



Mesa et a. 



torn [5] 



mm \3 



Macintosh et al. 



Macintosh et all (120061) [6] 
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Model 


EfTecti\'e Surfac-e 


S])()t Grou]) 1 


S])ot Grou]) 2 


Tyi)e 


Teui]) 


Ty])e 


Teuii) 


''iy])e 


Touii) 


Al 


Clear 


Ti 


Cloudy 


Ti 


Cloudy 


Ti 


A2 


Clear 


Ti 


Cloudy 


T2 


Cloudy 


T3 


Bl 


Cloudy 


Ti 


Clear 


Ti 


Clear 


Ti 


B2 


Cloudy 


Ti 


Clear 


T2 


Clear 


T3 


C 


Clear 


Ti 


Clear 


T2 


Clear 


T3 


D 


Cloudy 


Ti 


Cloudy 


T2 


Cloudy 


T;5 



Table 3: Six different realizations of the three distinct surface types (Effective Surface, Spot 
Group 1, Spot Group 2) covering the atmosphere of a giant planet: Al and A2) cloudy spots 
on a clear surface; Bl and B2) clear spots on a cloudy surface; C) clear surface with cold 
and hot clear spots; and D) cloudy surface with cold and hot cloudy spots. Ti, T2 and T3 
represent three different temperatures. 



